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This paper is concerned with the first and second passage times of a 
stationary random process, I(t, a), consisting of a sinusoidal signal of 
amplitude (#a) } plus stationary Gaussian noise with a finite expected 
number of zeros per unit time. This type of random process is present 
at the output of the IF amplifier of a radio or radar receiver during the 
reception of a sinusoidal signal immersed in Gaussian noise. Approximate 
integral equations are developed whose solutions yield approximate proba- 
ability densities concerning the first and second passage times of I(t, a). 
The resulting probability functions are presented in graphs for the case 
when the frequency of the sine wave is located in the center of a band of noise. 
Related results concerning the approximate distribution function of the 
absolute minimum or absolute maximum of I(t, a) in the closed interval 
[0, t] are also presented. 

I. INTRODUCTION 

Exact, explicit, results concerning the first passage times of a Markov 
or "Markov-like" random process have been given by many authors. 1-7 
But very little is known about the first passage times of a random 
process consisting of a sinusoidal signal plus stationary gaussian noise. 
This random process is of interest because it serves as a realistic model 
for the output of the IF amplifier of a typical radio or radar receiver 
during the reception of a sinusoidal signal immersed in Gaussian noise. 

Let I(t, a) denote the stationary random process consisting of a 
sinusoidal signal of amplitude (2a) * and angular frequency q plus 
stationary gaussian noise, I N (t), of zero mean and unit variance. Thus, 

I(t, a) = (2a)* cos (qt + 6 ) + I N (t). (1) 

do denotes a random phase angle which is distributed uniformly in 
the interval ( — ir, 7r). a denotes the signal-to-noise power ratio. 
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The first and second passage times of I(t, a) are indicated in Fig. 1 
and are defined as: 

(j) t + represents the time I(t, a) takes in going from an upcrossing 
of the level 7i to the first crossing of the level I 9 < I x . 

(u) t~ represents the time I(t, a) takes in going from a downcrossing 
of the level I x to the first crossing of the level I % < I t . 

(iii) t\ represents the time I(t, a) takes in going from an upcrossing 
of the level 7 X to the second crossing of the level / a < 1%. 

(iv) t~ x represents the time I(t, a) takes in going from a downcrossing 
of the level Ii to the second crossing of the level I 2 < It . 

For fixed I t , I 2 < I\ , and a we denote the probability densities of 
t + , r~, t\, and t" by W + (t, I lt I 2 , a), W~(t, I, , 7 2 , a), W\(t, I x , U , a), 
and W\(t, I x , I 2 , a), respectively. These four probability densities 
arise in many branches of science and technology. 

Because the random process I(t, a) is symmetrical about its mean 
value of zero, we need only discuss the case when I x > I 2 as is indicated 
in Fig. 1. The case when I x < I 2 can always be converted to the case 
under discussion by considering the random process — I(t, a). 

For the general process lit, a), exact, explicit expressions for the four 
probability densities are unknown. However, as already mentioned, 
exact, explicit results concerning the first passage times of I(t, a) are 
known for the special cases when a = and I N (t) is a Markov or Markov- 
like random process. 

The purpose of this paper is to present theoretical approximations for 
the four probability densities of the first and second passage times of 
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Fig. 1 — t + and r~ are the first passage times of I(t, a) = (2a) 1 ' 2 cos (qt + Bo) 
+Is(t) defined by the levels h and h ■ Similarly, n* and n" are the second 
passage times. 
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I(t, a) for the cases when a > and I N (t) is a stationary, gaussian process 
having a finite expected number of zeros per unit time. 

II. AUXILIARY PROBABILITY FUNCTIONS 

Using a notation consistent with Ref. 8, we define the following 
auxiliary probability functions concerning the stationary random 
process I(t, a): 

(0 PZ~(t, 1 1 , 1 2 , a)d,T, the conditional probability that a downward 
crossing of the level 7 2 occurs between t + r and t + t + dr given an 
upward crossing of the level 7, at t. 

(ii) PJ + (t, 7, , 7 2 , a)dr, the conditional probability that an upward 
crossing of the level 7 2 occurs between t -f- r and t + t + dr given a 
downward crossing of the level I t at t. 

(in) PJ~(t, Ix,I i} a)dr, the conditional probability that a downward 
crossing of the level I 2 occurs between t + t and t + r + dr given a 
downward crossing of the level /\ at t. 

(iv) PZ + (r, 1 1 , 1 2 , a)dr, the conditional probability that an upward 
crossing of the level I 3 occurs between t + r and t + t -\- dr given an 
upward crossing of the level /, at t. 

These auxiliary probability functions were given in Ref. S for the case 
when 7i = 7 a . Here, we need to merely generalize to the case when 
i\ 5^ I 2 . The reader should see Rice's work for the definition of all 
notation which is not defined in this paper. 4 When a > and 7, ^ 7 a , 
Rice's equation 38 generalizes to: 

P*~{t, 7i , 7 2 , a) 

= -pfcriW.]-' f dd f d/{ f dliWl V (h , I{ , U , h) (2) 

''-T ''O •'-00 

where 

N ti = Rice's equation 2.7 of Ref. 9 for the expected number of up- 
crossings (or downcrossings) of the level 7 X per second 

= Q8)» exp (-/?/2) A (-l)"(2n)l fa 
2tt ^o 2 n (n!) 3 \2 

•A(-i|;^(-|;» + i;-^) 

iF, = the confluent hypergeometric function 
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•exp l~ [M 22 (/{ 2 + J 2 2 ) + 2M i2 rJ[I' i + 2DJ[ + 2E 2 I' 2 + F 9 ]| 

2> a - if mCTi - Q cos 0) + Af l8 [Q cos (gr + 6) - I 2 ] 

+ iW 22 Qg sin + M 23 Qq sin (gr + 0) 

E a = M 12 [Q cos (gr + 0) - 7 8 ] + M 13 (/i - Q cos 5) 

+ M 22 Qg sin (gr + 0) + M 23 Qg sin 

F a = M n [Il + H- 2Q[7, cos + h cos (gr + 0)] 
+ Q 2 [cos 2 + cos 2 (gr + 0)]} 

+ 2M 12 Qq{ [I, - Q cos 0] sin + [Q cos (qr + 0) - I 2 ] sin for + 0) } 
+ 2M 13 Qg{[/i - Q cos 0] sin (gr + 0) + [Q cos (gr + 0) - I 2 ] sin 0} 
+ 2Af M [7, - Q cos 0][/ 2 - Q cos (qr + 0)] 

+ M 22 (Qg) 2 [sin 2 + sin 2 (gr + 0)] + 2M 23 (Qg) 2 sin 0sin (gr + 0). 
The M's are given in Ref. 4, Appendix I with 

m {j) = r W(f) cos 27r/r d/, (3) 

Jo 

where TF(/) = one-sided power spectral density of I N (t). Also, 
= — m"(0). The primes denote differentiations. 

Equation 2 can be put in the form: 

P» _ (r, h,h,a) = [4 7 r 2 AT /l ]- 1 M 22 (l - m 2 )' § 

■T exp (-G a /2M)J(n , h , fca) <*0 

where 

/(r, , *. , w - 27r(1 L ^ £ dx £ «? - w - * 3 k 



(4) 



a; 2 + ?/ 2 — 2risy 

1 ~ ~ 2(1 - r 2 ) 

fc, = M 22 [l - r\y\D 2 - 
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k 3 = -M" 22 [l - r^-\E 2 - ^[^jj^-J 

G 2 = M- 2 \[l - r^-\2nD 2 E 2 - D\ - E\] + F 2 . 

P 2 + (t, I x , 1 2 , a) is obtained from equation 2 by changing the signs 
of the co 's in the limits of integration. We find that P 2 + (t, I X1 J a , a) is 
equal to the right side of equation 4 with h 3 , k 3 replaced by — h 3 , —k 3 . 

P 2 ~(t, 7\ , 1 2 , a) is obtained from equation 2 by changing the upper 
limit of integration of I[ to — co . We find that P 2 ~(t, I x , I 2 , a) is equal 
to the right side of equation 4 with the function J(r x , h 3 , k 3 ) replaced 
by the function J^r, , h 3 , A; 3 ), where 

Ji(ri , h 3 , k 3 ) = jt-t: -2^ / dx I dy{x- h 3 )(y - k 3 )e'. (5) 

2tt(1 — fly J h , J k , 

P 2 + (r, Ij , 1 2 , a) is obtained from equation 2 by changing the lower 
limit of integration of I' 2 to + °° . We find that P 2 *(t, ?i , h , a) is equal 
to the right side of equation 4 with the function J(i\ , h 3 , k 3 ) replaced 
by the function «/i(r x , — h 3 , —k 3 ). 

The functions J(i\ , h 3 , k 3 ) and Jxfa , h 3 , k 3 ) are expressed in terms 
of Karl Pearson's well-known tabulated function (d/N) in Ref. 10. 

A considerable simplification occurs when a = 0. For this case we find 

P 2 + "(r, I, , J 2 , 0) - M 22 /T*(l - m 2 )" 5 

f"/ 2 (7 2 + 7 2 - 2^.^/2) 1 ,, , , x , A> 
' exp L 2 ~ 2(1 - m 2 ) J /(r ' ' *• ' k3) (6) 

where 

, —m'jmli — I 2 ) 

h3 " [M 22 (l - m 2 )]* 

, m'(/ i — m7 2 ) 

ft3 [Af M (l - m 2 )]*" 

Equation 6 reduces to Rice's 4 equation 47 when 7\ = 7 2 = 7. 

III. APPROXIMATE RESULTS VIA INTEGRAL EQUATIONS 

3.1 Probability Densities 

Let us assume that each of the random variables t + , t", t\ , and t\ 
(see Fig. 1), is statistically independent of the sum of the following 
(2AT + 2) axis-crossing intervals at level 7 2 when N = 0, 1, 2, • • • . Under 
this "quasi-independence" assumption, approximate theoretical results 
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for the probability densities of t + , t~, t\ , and t\ , namely W + (t, 7, , 7 2 , a), 
W-(t, 7, , h , a), W\{r, 7, , 7 2 , a), and W^t, 7, , I, , a) are given by the 
following integral equations: 

P + 2~(r, h,h,a) = W\r, h,h, a) 

+ W+(t, 7, ,h ,a)* P;~(t, h , h , a) (7) 
PVir, /, ,h,a) = W'(t, h , I 2 , a) 

+ W~(t, 7, ,h , a) * PT(t, 7 2 , / 2 , a) (8) 
P 2 ++ (r, /, , 7 2 , a) = TTftT, 7x , 7 2 , a) 

+ TF+(r, 7, , 7 2 , a) * P 2 + + (r, 7 2 , 7 2 , a) (9) 
P 2 " + (r, 7, , 7 2 , a) = T71(r, 7 X , I 8 , a) 

+ TTKt, A , U , a) * P 2 ++ (r, 7 2 , 7 2 , a). (10) 
The P 2 's are the auxiliary probability functions presented in Section 
II, and the symbol * denotes the convolution operator, that is, 

/ * Q m f 1(t)g(r - t) dt. (ID 

J — 00 

From symmetry we have 

PV(t, h,h,a) = P;-(t, h,h, a). (12) 

Integral equations 7 through 10 are analogous to McFadden's 
equation 47 and Rice's equation 84 in Refs. 11 and 4, respectively. 

Let us define two additional probability densities denned by the 
first and second passage times of I{t, a) between the levels h and 7 2 : 

(i) W{t, Ii,h, o)dr, the conditional probability that the first cross- 
ing of the level 7 2 occurs between t + t and t + t + dr given a crossing 
of the level h > 7 2 at t. 

(u) TPi(t, h, h, a)dr, the conditional probability that the second 
crossing of the level 7 2 occurs between t + r and t + t + dr given a 
crossing of the level 7 X > 7 2 at t. 
Clearly, we have 

W(r, I lt I 2 ,a) = \W\r, 7, , 7 2 , a) + £W~(r, I x , U , a) (13) 

W l {r, /,,/,, a) = jmCr, 7, , J a , a) + JlFfr, 7, , 7 2 , a). (14) 

3.2 A&soZwte Minimum In a Closed Interval 

Let us define the following probability functions concerning the 
absolute minimum of I{t, a) in a closed interval [0, t] : 
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F\t, h,h, a) 

s Pr { min /(/, a) > 7 2 | 7(0, o) -/,>/, , I'(0, a) > Oj (15) 

Oilir 

F-(t, h,h, a) 

= Pr { min I(t, a) > I 2 \ 7(0, a) = U > h , 7'(0, a) < 0| (16) 

0S*St 

F(r, 7, , I 2 , a) ^ Pr { min 7(*, a) > 7 2 | 7(0, a) = h> h\ (17) 

OSISr 

where Pr{ • } denotes the probability of the event inside the brace. 
Clearly, we have 



F + (r, h,h,a) = JT W\r, /,,/,, a) dr 

- 1 - f TF + (r, 7, , 7 2 , a) dr 

Jo 

F'(r, /,,/,, a) = JT W-(t, 7, , 7 2 , a) dr 

= 1 - f FT(r, A , U , a) dr 
F(r, 7, , 7 2 , a) = f~ W(t, 7, , 7 2 , a) dr 

- 1 - [ W(t, 7l , 7 2 , a) dr. 



(18) 



(19) 



(20) 



Because we are discussing only the case when 7i > 7 2 as is indicated 
in Fig. 1, we discuss only the probability functions concerning the 
absolute minimum of I(t, a) in a closed interval [0, t]. The corre- 
sponding probability functions concerning the absolute maximum of 
I(t, a) in a closed interval [0, t] are associated with the case when 
7i < I 2 , and they can be obtained from symmetry by considering 
the random process —I(t, a) . 

IV. RESULTS FOR SINUSOIDAL SIGNAL CENTERED IN LOW-PASS NOISE 

For purposes of computation we set the angular frequency, q, of the 
sinusoidal signal in the center of a band of gaussian noise with an 
ideal low-pass power spectral density of cutoff frequency / . Thus, 

q = 7r/ (21) 
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and 

W (f) - j^ 1 ° * ' f U (22) 

tOj othenvise. 

Accordingly, from equation 3, 

, . si n 2tj t /qo\ 

mW - -5ST' < 23 > 

From equation 23 we see that it is convenient to define normalized 
time as u = 2tt/ t. All our results are plotted with respect to normal- 
ized time wo. 

4.1 Experimental Verification of Auxiliary Probability Functions 

The auxiliary probability functions for the case when Ii = 7 a = are 
useful for studying the zero-crossing intervals, the axis-crossing intervals 
defined by the level I, = I 2 = 0, of I(t, a). Figures 2, 3, and 4 present 
P*~ , P;+ , P* + , and P 2 ~" for the case when Z, = h = and a = 0, 1, 
and 4. These results were computed by using Simpson's rule. The results 
compare satisfactorily with the initial behavior of the experimental 
probability densities presented in Figs. 34, 35, 42, 43, 46, and 47 of Ref. 
12. Notice that the experimental probability densities pertain to the case 
when the power spectral density of the noise is 

W (f) = 77TH (24) 

l + (£) 

rather than the power spectral density defined by equation 22. 

4.2 Results When a = and a = 1 

Figures 5 through 13 present the results when a = 0, signal absent, 
and h = 1, h = 0; h = 1, h = -1; and h = 0, I 2 = -1. The P 2 's 
were computed by using Simpson's rule, the integral equations defin- 
ing the W's were solved numerically by using the trapezoidal rule, 
and the F's were computed by using Simpson's rule. 

Similarly, Figs. 14 through 22 present the results when a = 1, signal 
present, and h = 1, h = 0; h = 1, I 2 — — 1; and I x = 0, I 2 = — 1. 

As an example of the interpretation of these results, we see from 
Fig. 16 that the median time, r«, for the random process I(t, 1) to 
go from the level I x = 1 to the level I 2 = for the first time is given by 

Wo = 2%f T m = 7T (25) 
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a =o 

1,-0 m <T,= Ji^ 




Fig. 2 (also Figs. 3 and 4) — Plots of the probability functions P a + ~(uc, h, la, 
a), Pi~*(iu>, h, h, a), Pa*{uo, h, h, a), and P a ~(.Uo, h, h, a) associated with 
the crossings of the levels h and I 3 by a stationary random process consisting of 
a sinusoidal signal of frequency / /2 plus stationary gaussian noise with autocor- 
relation function m{r). a denotes the signal-to-noise power ratio. 
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Fig. 3 — (See Fig. 2.) 
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Fig. 4— (See Fig. 2.) 



I. = i m(T)= 



sin u 




,A 



Fig 5 (also Figs. 8, 11, 14, 17, and 20) — Plots of the probability functions 
P.-(uo, h, h, a), W-(uo, h,h, a), P«*-(ito, h, h, a), W+(.iu>, h, I*, a), and 
W(uo, h, h, a) associated with the crossings of the levels h and I» bv a sta- 
tionary random process consisting of a sinusoidal signal of frequency / /2 plus 
stationary gaussian noise with autocorrelation function m{r). a denotes the 
signal-to-noise power ratio. 
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Fig. 6 (also Figa. 9, 12, 15, 18, and 21)— Plots of the probability functions 
Pr*(ifa, h, I-., a), TTr(uo, h, h, a), JY + (uo, h, h, a), TPi + (uo, It, h, a) and 
Wiiuo , h, h, a) associated with the crossings of the levels h and h by a sta- 
tionary random process consisting of a sinusoidal signal of frequency / n /2 plus 
stationary gaussian noise with autocorrelation function to(t). a denotes the 
signal-to-noise power ratio. 




U =27Tf T 

Fig. 7 (also Figs. 10, 13, 16, 19, and 22) — Plots of the probability functions 
F*(uo, h, h, a), F~(no, h, h, a), and P(uo, h , It, a) associated with the 
crossings of the levels h and h by a stationary random process consisting of a 
sinusoidal signal of frequency / /2 plus stationary gaussian noise with auto- 
correlation function mir). a denotes the signal-to-noise power ratio. 
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Fig. 8— (See Fig. 5.) 
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Fig. 9— (See Fig. 6.) 
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Fig. 10 — (Sec Fig. 7.) 
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Fig. 11 — (See Fig. 5.) 
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Fig. 12— (See Fig. 6.) 
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Fig. 13 — (See Fig. 7.) 
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Fig. 15— (See Fig. 6.) 
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U =21Tf T 

Fig. 16— (See Fig. 7.) 
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Fig. 17 — (See Fig. 5.) 
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Fig. 18— (See Fig. 6.) 
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Fig. 19— (See Fig. 7.) 
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Fig. 20 — (See Fig. 5.) 
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Fig. 21 — (See Fig. 6.) 
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Fig. 22 — (See Fig. 7.) 



or 



r m = (2/c)- 1 . (26) 

r m also represents the median time during which the random process 
I(t, 1) remains continually above the level I 2 = when it starts at the 
level 7i = 1. From symmetry, t„, also represents the median time dur- 
ing which the random process I(t, 1) remains continually below the 
level 7o = when it starts at the level Ii = — 1. 



V. CONCLUSIONS 

The exact auxiliary probability functions can be used in approxi- 
mate integral equations in order to deduce approximate probability 
densities of the first and second passage times of sine wave plus sta- 
tionary, gaussian noise with a finite expected number of zeros per 
unit time. These approximate probability densities can be used to 
deduce the approximate median times associated with the first and 
second passage times. Also, the approximate probability densities of 
the first passage times can be used to deduce approximate distribu- 
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tion functions for the absolute minimum or absolute maximum of sine 
wave plus noise in a closed interval. 

The corresponding exact results are not yet known. 
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